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Contributions of different parts of spin-spin interactions to quantum correlations in a 

spin ring model in an external magnetic field 
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We study quantum correlations in a bipartite heteronuclear (N — 1) x 1 system in an external 
magnetic field. The system consists of a spin ring with an arbitrary number IV — 1 of spins on 
the ring and one spin in its center. The spins on the ring are connected by secular dipole-dipole 
interactions and interact with the central spin through the Heisenberg zz-interaction. We show that 
the quantum discord, describing quantum correlations between the ring and the central spin, can 
be obtained analytically for this model in the high temperature approximation. The model allows 
us to find contributions of different parts of the spin-spin interactions to quantum correlations. We 
also investigate the evolution of quantum and classical correlations at different numbers of spins. 

PACS numbers: 03.67.-a, 76.60.-k 


INTRODUCTION 


Quantum correlations in many-qubit systems are responsible for the effective work of quantum devices (in particular, 
quantum computers) and give them significant advantages over their classical counterparts [l|. In order to create 
quantum devices we have to study quantum correlations and to control them. Spin systems are very suitable for the 
study of quantum correlations due to the developed theoretical methods Q and the experimental methods of nuclear 
magnetic resonance (NMR) 0] for the excitation, detection and coherent control of the individual spin states. 

Until recently it has been accepted that entanglement is responsible for quantum correlations, and quantum devices 
can be created only on the basis of using materials with entangled states Q. However, it has turned out that 
quantum algorithms ;J| which significantly outperform the classical counterparts can work using mixed separable 
(non-entangled) states. Furthermore, it turned out that quantum non-locality can be observed in systems without 
entanglement ’7J. From this we can conclude that entanglement describes only a part of quantum correlations but not 
all of them. According to the current understanding [5], total (quantum and classical) correlations in a system are 
defined by the mutual information. The problem is how to separate the classical correlations from the quantum ones. 
This problem was solved independently by Henderson and Vedral [S] as well as Ollivier and Zurek j§]. The classical 
correlations in a two-partite system are determined by a complete set of projective measurements carried out only 
over one of the subsystems [§j. Then a measure of quantum correlations (the quantum discord) is determined as the 
difference between the mutual information and its classical part, maximized over all possible projective measurements 
fil, [ 9 ]. The quantum discord is determined completely by quantum properties of the system and equals zero for 
classical systems. 

Computing the quantum discord is a rather tedious problem because we must optimize the quantum conditional 
entropy over all possible complete sets of measurements in any subsystem of the system under study [§, [ij[ . It 
was shown fiol ] that computing quantum discord is NP-complete. The already developed methods [11—13] permit 
computing the discord only in two-qubit systems [14—17] and in the simplest three-qubit systems SGI. At the 
same time, it is very important to calculate the quantum discord in many-qubit systems. For example, the many-qubit 
quantum discord is significant for solving problems of NMR quantum information processing Q. It is well known that 
NMR is used to obtain an experimental implementation of quantum algorithms in few-qubit systems. The quantum 
entanglement vanishes at room temperatures in highly mixed states that occur in liquid-phase NMR experiments 
[2(] ]. This raised doubts that NMR can be used to demonstrate the advantages of quantum algorithms over their 
classical counterparts in multiqubit systems Q. The quantum discord does not vanish 21] in mixed states in NMR 
experiments in contrast to the entanglement and indicates that quantum correlations exist in many-qubit systems. 
Experimental studies [21] using NMR-tomography methods ![J] support this conclusion. The development of methods 
for computing the quantum discord in many-qubit systems is therefore very relevant. 

There are a lot of different NMR methods in solids which can be used for experimental investigations of quantum 
correlations. First of all, notice the methods of spin echo and free induction decay [ 2 , Q] ■ The time evolution of 
signals of spin echo and free induction decay in the system of an electron (or a hetero-nucleus) surrounded by a cloud 
of other nuclei can give very important information about quantum discord and its evolution [22|. We notice also 
multiple quantum NMR methods 23] which allow us to obtain the time and temperature dependencies of quantum 
correlations 124 1. 
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FIG. 1: Many-spin systems for the investigation of the quantum discord. Subsystem A is a ring of spins, coupled by the DDI. 
The impurity spin (in the center of the ring), connected by the ^-interactions with the ring spins, represents subsystem B. 


We have suggested 19j, 251 a spin model for the calculation of the quantum discord. In our model, a linear chain 
of spins, coupled by the dipole-dipole interactions (DDI) in a strong external magnetic field, is connected with an 
impurity spin by the Heisenberg ^-interaction. This model can be used for the analysis of NMR experiments in solids, 
both at low and high temperatures. The developed approach fl9l . [25 ] allows us to study quantum correlations and 
their unitary evolution in the conditions of NMR experiments. Unfortunately, in that model the quantum discord can 
only be calculated numerically with the random mutation algorithm 26}. The analytical calculation of the quantum 
discord for the model 0 is possible only for three-qubit systems [25]. Later we understood |2?]| that it is possible to 
find the quantum discord analytically for many-qubit systems if one considers a spin ring instead of a spin chain used 
in IS, |25(. We note that this model can be applied not only to heteronuclear systems, but also to electron-nuclear 
ones. If we consider a central electron spin surrounded by a ringof nuclear spins then the model is close to the model 
of deterministic quantum calculations with one qubit (DQC1) [6]. 

It is necessary to emphasize that the multiqubit spin ring model can be investigated by standard tools [11—13] 
for two-part binary systems. However, we simplify those methods with an optimization of the quantum conditional 
entropy over the parameters of the projectors. As a result, this approach allows us to investigate contributions of 
different interactions to quantum correlations. 

The main goal of the article is the investigation of the quantum correlations in a spin ring system and the deter¬ 
mination of the contributions of different parts of spin-spin interactions to quantum correlations at different relation 
between the Larmor frequencies of the ring spins and the impurity one. 

The paper is organized as follows. A spin ring system, coupled with a central spin, in the external magnetic field 
is presented in Sect. 2. The optimization of the quantum conditional entropy in the high temperature approximation 
is considered in Sect. 3. Calculations of the quantum discord are given in Sect. 4 at different relations between the 
Larmor frequencies of the spin subsystems. A comparison of quantum and classical correlations in the course of the 
system evolution is given in Sect. 5. We briefly summarize our results in Sect. 6. 


A BIPARTITE MANY-QUBIT SYSTEM IN AN EXTERNAL MAGNETIC FIELD 

We consider a ring of nuclear spins (I = 1/2), coupled by the DDI in a strong external magnetic field. We will call 
it subsystem A. The distances between spins of subsystem A are not necessarily the same. Subsystem A interacts 
with an impurity spin S, which is in the center of the ring (subsystem B , Fig. 1). The strong external magnetic field 
is perpendicular to the plane of Fig. 1. 

We will consider the problem in the high temperature approximation [2]. Thes approximation is widely used in 
different problems of spin dynamics and magnetic resonance [2]. For example, the high temperature approximation 
in systems of spins coupled by DDI is valid [2. 28] at 


/3co a < 1 , fivs < 1 , ( 1 ) 

where [3 is proportional to the inverse temperature |2] of the system; uja and ojb are the Larmor frequencies. 

The proof is given in Ref. [28| and it is based an the fact that the DDI is inversely proportional to the cube of the 
distances between spins. However, it is not so in the considered case. For example, N — 1 spins interact with the 
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central spin and all coupling constants are equal to each other. The dependence of the DDI coupling constants of the 
ring spins on the distances between spins is very weak and the conditions © are insufficiently. It is evident that the 
conditions should be modified. We expect that the high temperature approximation holds when 


(N — 1)/3uia < 1) {N — 1)/3wb < 1. 


( 2 ) 


The conditions of Eq. © are valid for small spin clusters. Decoherence (in particular, external noise) is determined 
by the relaxation times Tj and T% @ which are sufficiently long for such clusters in NMR experiments [29, 30]. Thus, 
we can neglect noise in our spin systems in the calculation of quantum correlations. 

In the initial moment of time, the system is in the thermodynamic equilibrium state. In the high temperature 
approximation [ 2 ], the density matrix p(0 ) is 


p(0) — (1 + P^aIz + /3ubSz), 


(3) 


where z points in the direction of the external magnetic field; N is the total number of the spins (N — 1 is the 
number of the ring spins); I z = JT I z i, and I zz , S z are the z-projections of the i-th spin on the ring and the impurity 
spin. Applying resonance (for spins / and S) 90°-radiofrequency pulses of the magnetic field around the axis y in the 
rotating reference frames (for spins /*, i = 1,2,..., IV — 1 and S ) [2j, we create conditions for the emergence of quantum 
correlations during the evolution of the spin system. Formally, such initiating pulses lead to the replacement of the 
operators Ii Z , S z with Ii X and S x , where Ii a and S a (a = x, y, z) are projections of the i-th ring spin and the impurity 
spin on the axis a. The subsequent evolution of the system is described by the Hamiltonian H = Hd z + H zz . Here 
the Hamiltonian Hd z describes the secular DDI in the strong external magnetic field in the spin ring (subsystem A). 
This secular part of the DDI can be written as [ 2 ] 

Hdz = ^2 dij (siizljz - IJj) > ( 4 ) 

i<j 


where dij is the DDI coupling constant of ring spins i, j: and Iilj = Ii X Ij X + hyljy + hzljz- The Hamiltonian H zz 
characterizes the interaction between subsystems A, and B and equals 

H zz = Y j ahzS z , (5) 

i 


where g is the coupling constant of the ring spins with the impurity spin, which has a different gyromagnetic ratio from 
the ring spins. We took into account that the coupling constants of the zz-interactions in the considered case (see Fig. 
1) are the same for all ring spins (see Eq. ©). It means that the Hamiltonians Hd z and H zz commute at an arbitrary 
number of the ring spins. This property holds in the model with a linear spin chain instead of a spin ring only for a 
chain with two spins fl9|. Thus, subsystem A is only subject to the dipole evolution in the considered case. Since a 
unitary local transformation does not change the discord, we can ignore the dipolar evolution of subsystem A in our 
calculations. Taking into account that the evolution of the system is determined by the Heisenberg zz-interactions 
only, one can write the density matrix of the system as 

p(r) = e~ lHzzt p(0)e lH ‘* t = 277 {l + /3u} A [Ix cos( r) + 2 I y S z sin(r)] 

1 1 ( 6 ) 
+/3lu b [S x cos(2 tI z ) + S y sin(2r/,)]|, 


where t = gt/2 is the dimensionless time. 

The reduced density matrices pa(t), Pb(t) are necessary for the calculation of the quantum discord. They can be 
represented as 


Pa{t) 


1 

2^-1 


[1 


+ /3u a I x cos(t)], 


(7) 


Pb(t) = 5 jl + 27 v=t/ 3 wb S' x Tr[cos(2r I z )] + S' y Tr[sin(2r/,)]J j 


— 2 [1 + f3u B S c 


cos N 1 (r) 


We used in Eq.© the relationships 

Tr[cos(2 tI z )} = 2 JV “ 1 cos n ~\t), Tr[sin(2r/ 2 )] = 0. 


( 8 ) 


( 9 ) 






These relationships © are proved in Appendix A. 

We will also use the expression for the entropy S(p) of the system 

S (P ) = n ~s^ 2[( n - 1 )P 2 “a + P 2 “ 2 b\ 

and the expressions for the entropies S(pa) and S(pb) of subsystems A and B are 
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( 10 ) 


S(p A ) = N-1- -^y/3 2 ^ cos 2 (t), 


(ID 


S(pb) = 1 - cosHN 


( 12 ) 


THE OPTIMIZATION OF THE QUANTUM CONDITIONAL ENTROPY IN THE HIGH 

TEMPERATURE APPROXIMATION 


According to the standard approach 11, 12| for the optimization of the quantum conditional entropy one performs 
a total set of projective measurements over one-qubit subsystem B, { B k = VA k V', k = 0,1}, where the matrix 
V € SU(2) and II k (k = 0,1) are projectors. For our aims, it is better to write the projectors through the spin 
operators 


n 0 - 2 A Tt X S X 4“ TlySy A TI Z S Z1 III - 2 B X S X BySy n z S z , 


(13) 


where n x , n y , n z are the coordinates of the unit vector in the spin space (n 2 +n 2 + n 2 = 1). It is evident that A 2 = life 
(k = 0,1). Using the projectors (THU) we can replace an arbitrary unitary matrix V € SU{2). At the same time, such 
projectors are more convenient for finding contributions of different parts of the DDI to quantum correlations. 

The following relation is very useful at further calculations 


(~l) fc 

n fe s a n fe — — n a A k: 


a = x,y, z; k = 0,1. 


(14) 


The density matrix p(t) of Eq. (jGj can be transformed after performing measurements with Eq. m as 

n k p(t)n k = 27 V \l A +/3wa[4cos(t) + {-l) k n z Iy sin(r)] 

+{—l) k \fiu)B [n x cos(2r I z ) + n y sin(2r/ 2 )] | ®II fc , k = 0,1, 


(15) 


where I a is a unit operator in the spin space of subsystem A. Thus, one can find that the whole system is described 
by the ensemble of the states {p k ,p k }(k = 0,1) after the measurements, where 


and the matrices po , p\ are 


p 0 = Tr{n 0 /9(t)n 0 } = \ + jn x /3u} B cos N 1 (t), 
Pi = Trjnip^jlli} = \- \n x f}u> B cos^ -1 ^), 


Po = 2 ^ {l + P^a[I x cos(t) + n z I v sin(r)] 
+ 1 / 3 u) B [n x cos(2 tI z ) + n y sin(2r/ 2 )]|, 
pi = 2 ^{l +/3wn[4cos(r) - n z I y sin(r)] 
— \P u B [n x cos(2 tI z ) + n y sin(2r/ 2 )]|. 


(16) 


(17) 


The results of Appendix A are used again at calculations of po and p± in EqS; (fTBl) . Then the conditional quantum 
entropy S con d after the measurements over subsystem B can be written llj,jl2| as 


* 5 cond — PO^(po') + PlS'^Pl'): 


(18) 
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where S(pk ) = —Tr[p^ log 2 pk]- It is convenient to introduce the parameters u and v, describing the high temperature 
approximation (see Eq. ©) 


U = - 1)/?WB < 1, V = ^(N - 1)0WA < 1. (19) 

The calculation of the quantum conditional entropy Cl with equations CSbCIlb d up to the terms of the orders 
u 2 and v 2 leads to the formula 


Scond — 2 In 2(iV—l) 2 { 7 


2 1+cos jv ~ 1 (2t)-2cos 2( - w ~ i:> (t) 


— (N — l)v 2 sin 2 (r) 


+nt 


2 1 cos 2 —— (N — l)v 2 sin 2 (r) + a(u,u)j, 


( 20 ) 


where the function a(u, v) = (N — \)v 2 — 2 (N — l) 3 In 2 does not depend on n a (a = x, y, z). 


THE QUANTUM DISCORD AT DIFFERENT RELATIONS BETWEEN THE LARMOR FREQUENCIES 

OF THE SPIN SUBSYSTEMS 


In this section we demonstrate that quantum correlations depend on relations between the Larmor frequencies of 
the ring spins and the central one. First, we consider the case of v > u, i.e., the Larmor frequency of the ring spins 
exceeds those of the central spin. We show in Appendix B that the coefficients of equation (1201) at n x and n y are 
positive in this case. It means that the quantum conditional entropy of equation (1181) achieves the minimal value at 


n x = n y = 0, \n z \ = 1. 


( 21 ) 


In this case (see Eq. (1151) 1 only I y ^-interactions can give a contribution to quantum correlations due to the term 
which is proportional to n z I y . 

Using Eqs. cud, cn m and m, one can obtain an expression for the quantum discord D , which is the 
entropic measure of the quantum correlations 


D = 


21n2(./V — l) 2 


l-cos 2 ^- 1 ^) 


( 22 ) 


Naturally, the quantum discord here is determined by the lower Larmor frequency u of the impurity spin. We notice 
also that formula J22D is valid for such numbers of the ring spins when the conditions © of the high temperature 
approximation Q, l28l] hold. 

The dependence of the quantum discord on the dimensionless evolution time r at the different numbers of the spins 
with u/(N — 1) = \(3ub = 0.015 is shown in figure 2. 

The case u > v is very difficult for an analytical investigation. However, it is possible to perform the optimization 
of the quantum conditional entropy at some additional conditions. Let 


u 2 >(N-1)v 2 , 0 < t < 7t/4, TV — 1 > 2. 


(23) 


Then we prove in Appendix C that the quantum conditional entropy of equation (1181) achieves the minimal value at 


n x =n z = 0, 


Kl = i- 


(24) 


We can conclude again from equation (USD that -UiSy-interactions are responsible for quantum correlations at the 
conditions (1^51) . Using Eas. lflOl) . (fl2l) . (P2U1) and ©j}, one can calculate the corresponding quantum discord: 


D = 


2 In 2 (N 


jv-ip ~ 1 ) t;2 sin 2 (r) + ^ cos N 3 (2t) + 1 - 2 cos 2(w 1) ( T )]}- 


(25) 


It is worth to notice that I z S y - interactions are also responsible for quantum correlations not only at the conditions 
(121 but also when 


u 2 > 


N — 1 


3 1 - 3-Lv-!) 


2 , 7t/4 < t < arctan(-\/2), 


(26) 


Then again the quantum conditional entropy achieves its minimal value when Eq. (1241) holds, and the quantum 
discord is determined by formula (1251) . 
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FIG. 2: The dependence of the quantum discord on the evolution time at different numbers of the spins (N = 3,7,11) with 
u = 0.015. The minimal value of the conditional entropy is realized at n x = n y = 0 and \n z \ = 1. 



FIG. 3: The regions of the parameters 7 = uja/^b and the dimensionless evolution time r, where I y S z ~, I z S y -, I Z S X - spin-spin 
interactions are responsible for quantum correlations at IV = 5. 


Finally, we show in Appendix C that at the conditions 

u 2 > 2(TV — \)v 2 , N > 3 is odd, arctan(\/2) < r < 7 t/ 2 , (27) 

the quantum conditional entropy of equation (1181) achieves its minimal value at 

n y =n z = 0 , \n x \ = 1 . (28) 

We find with equation m that /siS^-interactions are responsible for quantum correlations in this case and the 
quantum discord is 


D = 


1 

2 In 2 (AT — 1) 


\)v 2 sin 2 (r) + ^ 1 — cos^ 1 (2r) j 


(29) 


In Fig. 3 we demonstrate the regions, where I y S z ~, I z S y - and /^^-interactions are responsible for quantum 
correlations at different 7 = uja/wb and the evolution time r for N = 5. 












7 




FIG. 4: A comparison of quantum and classical correlations in the course of the evolution of the spin system at N = 9; a) an 
initial period of the evolution 0 < r < 7t/4; b) the evolution on the time interval arctan (%/2) < r < tt/2. 


A COMPARISON OF QUANTUM AND CLASSICAL CORRELATIONS IN THE COURSE OF SPIN 

EVOLUTION 


Classical correlations can be obtained from the maximal information about subsystem A after measurements on 
subsystem B [Bj, [gj]. The total correlations I(p) are determined by 


I(p) = S(p A ) + S(p B ) ~ S(p), 


(30) 


where the entropies S(pa), S(pb), S(p) are determined by the expressions (fTUl) . ([TT1) . (fl2l) . Taking into account 
formulas ([23), ([23, (EH) for the quantum discord, one can obtain the expressions for classical correlations C at the 
different conditions for the minimal value of the quantum conditional entropy: 


C = 


1 


2 In 2(iV — 1) 


! sin 2 (r), for 


= 1 , 


^y 0 , 


(31) 


_ U 

~ 41n2(JV—l)* 
for \n x \ = 1, 


1 + cos^ 1 (2r) — 2 cos 2 ^ w 1 ^(r) , 
riy — ti z — 0, 


(32) 


C = 


41n2(AT- l) 2 L 


1 — cos^ 1 (2 t ) 


for 


= 1 , 


= n, = 0. 


(33) 


We can compare quantum and classical correlations in the course of the evolution at different numbers of spins. 
Such comparison is presented for the system, consisting of N = 9 spins, in Fig. 4. Fig. 4a demonstrates that classical 
correlations exceed quantum ones at the initial period of the evolution. However, the quantum correlations increase 
with time and exceed classical at 0.521 < r < 7r/4. The evolution of quantum and classical correlations on the time 
interval [arctan (y/2), 7r/2] is given in Fig.4b. Quantum correlations exceed classical ones on the main part of this 
interval. However, quantum correlations rapidly decrease and classical correlations increase at the end of the interval 
[arctan ( 1 / 2 ), tt/2]. 


CONCLUSIONS 

Our work is aimed at the analytical computation of the quantum discord in a system of interacting spins in a strong 
external magnetic field. The spin dynamics of this system can be investigated experimentally by NMR methods. 



















Our results demonstrate that there is a non-zero quantum discord at room temperature in the considered spin 
system. It is relevant for NMR in liquids where entanglement is absent at room temperature. Our results mean that 
there are quantum correlations in the considered system at room temperature and possibilities of NMR methods for 
creating quantum devices and performing quantum calculations are not exhausted. 

We modified the standard methods ll|, 121 for the calculation of the quantum discord in bipartite systems. Our 
idea is to use projectors for quantum measurements instead of arbitrary unitary transformations in order to optimize 
the quantum conditional entropy over parameters of those projectors. 

We showed that different parts of the spin-spin interactions are responsible for quantum correlations at different 
relations between the Larmor frequencies of the ring and impurity spins. It is very important that the model considered 
in this paper allows us to calculate the quantum discord analytically and we obtained the analytical formulas for the 
quantum discord at arbitrary relations between the Larmor frequencies of the ring spins and the impurity one. We 
have also investigated the time evolution of the quantum discord in the spin ring system with a central spin in the 
external magnetic field. The quantum and classical correlations were compared in the course of the spin system 
evolution. 

Our model is also likely to be useful in the investigation of quantum and classical correlations at low temperatures 
where the quantum discord is significantly larger than at room temperature. 

The work is supported by the Russian Foundation for Basic Research (Grants No. 13-03-00017, No. 13-03-12418 
and No. 15-07-07928) and the Program of the Presidium of RAS No. 32 “Electron spin resonance, spin-dependent 
electron effects and spin technologies”. 


Appendix A 

The ring contains TV — 1 spins and one can write that 

N-l 

Tr{cos(2r/ z )} = ^ Cn_ 1 cos 
fc=o 

where Cjy_i = (^(T 1 ) = k\(N~i-k)\ ■ R * s suitable to use the complex values: 

Ek=o C k N - r cos [2 t - k)] = Re {E^ C k N _} 

= Re J2k=o C^„ ie - 2lTfe } 

Applying the binomial theorem one obtains 

Re EfE 1 C n- ie~ 2irk J = Re { e ir ( JV - 1 )(l + e- 2 ^” 1 } 

= Re {(e lT + e~ lT ) N ~ x } = 2 JV_1 cos n ~ 1 (t). 

Now consider Tr{sin(2r/ Z )} and perform the rotation e lnIx of the operator sin(2r/ z ): 

Tr{sin(2 tI z )} = Tr{e i7r/ * sin(2 Tl z )e~ inI *} 

= —Tr{sin(2r I z )}. 

It means that Tr{sin(2r/ Z )} = 0. The relations (0 are proved. 

Appendix B 

We need to prove (see Eq. d^Il l that 



(A34) 


(A35) 


(A36) 


(A37) 


1 + cos^ 1 (2r) — 2cos 2 ( w 1 ^(r) < 2(TV — 1) sin 2 (r). 


(B38) 
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The left hand side of Ea. (IB38|) can be rewritten as 

1 + cos Ar_1 (2r) — 2 cos 2 ^ -1 ) (r) 

= 1 — cos 2 i Ar_1 )(r) + cos 7V_1 (2 t) — cos 2 ( Ar_1 )(r) 

= sin 2 (r) (l + cos 2 (r) + cos 4 (r) + ... + cos 2 ^ -2 ) (t)) 

+ ^cos(2r) — cos 2 (r)j ^cos Ar_2 (2r) +cos Ar_3 (2 t) cos 2 (r) 
+... + cos 2 ^- 2 )^)) < sin 2 (t){N - 1) + sui 2 (t)(A - 1) 

The statement follows from the inequality (IB39I) . 

Further we need to prove (see Eq. m that 

1 — cos iv_1 (2r) < 2 (N — 1) sin 2 (r). 


(B39) 


(B40) 


We can transform the left hand side of Eq. (IB40I) in the following way 

1 — cos Ar_1 (2r) = ^1 — cos(2t)^ ^1 + cos(2r) + cos 2 (2t) 

+... + cos Ar ~ 2 (2 r)j < 2 sin 2 (t)(1V — 1) 

The statement (IB40D follows from the inequality (1B41II . The inequalities (IB38I) and (IB40D lead to the conditions (12T1) 
for the minimal value of the conditional entropy of Eq. El- 


Appendix C 


We will show that the coefficient of n 2 in the curly brackets of Eq. m exceeds the one of n 2 at the conditions 
(El. One obtains from Eq. (EOl) that we need to prove that 


2 1—cos'™ 1 (2 r) / A T i\ 2 • 2/ \ v 2 1+cos™ 1 (2r)—2cos' 

tr- 2 —'—- — {N — ljirsin (r) > u z — - 1 —^2 - 

— (N — \)v 2 sin 2 (r). 

It is evident that Eq. (IC42I) can be transformed to 

1 > [1 — tan 2 (r)] Ar_1 . 


2(IV — 1) 


( T ) 


(C42) 


(C43) 


Eq. (IC43I) is valid at the conditions E51) . Now we prove that the coefficient of n 2 in the curly brackets of Eq. El is 
non-negative at the conditions E31) . In fact, we should prove that 


1 — cos^ 4 (2r) 


> sin 2 (r). 


Eq. (IC44D is equivalent to 


cos(2r) > cos iv i (2r) 


(C44) 


(C45) 


and Eq. (IC45I) holds due to the conditions El. One can conclude from Eq. El that the minimal value of the 
quantum conditional entropy is achieved at the parameters 


Next, we will show that the coefficient of ro 2 in the figure bracket of Eq. El exceeds the one of n 2 at the conditions 
El- In this case we can obtain that this statement is valid, if 

[1 — tan 2 (r)] 2 > 1. (C46) 

However, it is really valid at the conditions El- We show now that the coefficient of n 2 in the curly bracket of El 
is non-negative. We would like to prove that 


2 1 —1— cos^ 1 (2r) - 2 cos 2(w x )(2t) 2 . 2 

u 2 --—--—- - (N - l)v 2 snr(r) > 0. 


(C47) 
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Using the conditions m, one can obtain from Eq. (IC47I) that 

[cos 2 (t) — cos^“^(r)] + [cos Ar_1 (2r) — cos 2 ^ jv_1 ^(t)] > 0. (C48) 

The expression in the first square bracket is non-negative at N > 3. The second square bracket can be transformed 
as 


cos 


N -1 


(2 t) — cos 2 ^ N 1 ^(r) = cos 2(JV 1 ^(r) f 1 — tan 2 (r)^ 


N-l 


- 1 


> 0 . 


(C49) 


This term is positive because TV — 1 is even and arctan (V2) < r < n/2 due to the conditions (l27l) . In fact, we proved 
that the quantum conditional entropy of (1181) achieves its minimal value at the conditions (051) . 


* Electronic address: s.i.doronin@gmail.com 
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